Meter electronics and fluid quantification method for a fluid being transferred

ABSTRACT

Meter electronics ( 20 ) for quantifying a fluid being transferred is provided. The meter electronics ( 20 ) includes an interface ( 201 ) configured to communicate with a flowmeter assembly of a vibratory flowmeter and receive a vibrational response and a processing system ( 203 ) coupled to the interface ( 201 ). The processing system ( 203 ) is configured to measure a volume flow and a density for a predetermined time portion of the fluid transfer, determine if the fluid transfer is non-aerated during the predetermined time portion, if the predetermined time portion is non-aerated then add a volume-density product to an accumulated volume-density product and add the volume flow to an accumulated volume flow, and determine a non-aerated volume-weighted density for the fluid transfer by dividing the accumulated volume-density product by the accumulated volume flow.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to a vibratory flowmeter and method, and more particularly, to a meter electronics and fluid quantification method for a fluid being transferred.

2. Statement of the Problem

Vibrating conduit sensors, such as Coriolis mass flowmeters and vibrating densitometers, typically operate by detecting motion of a vibrating conduit that contains a flowing material. Properties associated with the material in the conduit, such as mass flow, density and the like, can be determined by processing measurement signals received from motion transducers associated with the conduit. The vibration modes of the vibrating material-filled system generally are affected by the combined mass, stiffness and damping characteristics of the containing conduit and the material contained therein.

A typical Coriolis mass flowmeter includes one or more conduits that are connected inline in a pipeline or other transport system and convey material, e.g., fluids, slurries, emulsions, and the like, in the system. Each conduit may be viewed as having a set of natural vibration modes, including for example, simple bending, torsional, radial, and coupled modes. In a typical Coriolis mass flow measurement application, a conduit is excited in one or more vibration modes as a material flows through the conduit, and motion of the conduit is measured at points spaced along the conduit. Excitation is typically provided by an actuator, e.g., an electromechanical device, such as a voice coil-type driver, that perturbs the conduit in a periodic fashion. Mass flow rate may be determined by measuring time delay or phase differences between motions at the transducer locations. Two such transducers (or pickoff sensors) are typically employed in order to measure a vibrational response of the flow conduit or conduits, and are typically located at positions upstream and downstream of the actuator. The two pickoff sensors are connected to electronic instrumentation. The instrumentation receives signals from the two pickoff sensors and processes the signals in order to derive a mass flow rate measurement, among other things. Vibratory flowmeters, including Coriolis mass flowmeters and densitometers, therefore employ one or more flow tubes that are vibrated in order to measure a fluid.

Bunkering refers to the practice of storing and transferring marine fuel oils, which have come to be known as bunker fuels. For ship fueling, large amounts of fuel may be temporarily stored in a barge or other container for the purpose of transferring fuel from shore to a ship. A bunker may be located on a dock or other port facility, or may be carried by a barge or other refueling vehicle. During bunkering, the fuel measurement usually comprises an empty-full-empty batching process, allowing gas to become entrained in the fuel.

Prior art fuel oil bunkering methods are based on volumetric tank measurements and a reference density typically obtained by laboratory sample. Look-up tables and reference density measurements are traditionally used to calculate the total mass of the delivered bunker fuel, in conjunction with a tank level measurement or dip tape measurement. Prior art measurement accuracy depends on many factors, including temperature, pressure, the presence or absence of entrained gas, dip tape measurement error or uncertainty, tank volume uncertainty, accuracy of conversion tables, human error, and how well the density sample represents the average batch density, for example. Further, heavy fuel oil (HFO) tends to stratify over time, so components may separate and have different densities, viscosities, etc.

While mass flow total is the most important measurement for a fuel transfer operation, the density and viscosity are also desired, as they are important indications of the quality or grade of the fuel oil. Generally, the higher the density and viscosity the lower the quality of fuel delivered. Knowledge of density and viscosity is essential for ensuring that the ideal fuel for a particular engine is used. In addition, a customer will want to determine how much of a fluid transfer is aerated. Such a determination can indicate whether the fuel is being intentionally aerated in order to increase the apparent volume of the fuel being delivered.

In the prior art, the current method for determining the bunker fuel's density and viscosity is to take a sample during the bunkering operation and send it to a laboratory for analysis. Unfortunately, laboratory analysis typically takes several days and the fuel will be on the ship and being used by that time. Further, fuel oil tends to stratify in tanks, leaving the lower grade fuel with the highest viscosity and density at the bottom of the tank and the higher quality fuel at the top. A single sample cannot capture the average density or viscosity of an entire batch.

ASPECTS OF THE INVENTION

In one aspect of the invention, meter electronics for quantifying a fluid being transferred comprises:

-   -   an interface configured to communicate with a flowmeter assembly         of a vibratory flowmeter and receive a vibrational response; and     -   a processing system coupled to the interface and configured to:         -   measure a volume flow ({dot over (ν)}_(i)) and a density             (ρ₁) for a predetermined time portion (t_(i)) of the fluid             transfer;         -   determine if the fluid transfer is non-aerated during the             predetermined time portion (t_(i));         -   if the predetermined time portion (t_(i)) is non-aerated,             then add a volume-density product ({dot over (ν)}_(i)ρ_(i))             to an accumulated volume-density product ({dot over             (ν)}ρ_(accum)) and add the volume flow ({dot over (ν)}_(i))             to an accumulated volume flow ({dot over (ν)}_(accum)); and         -   determine a non-aerated volume-weighted density             (ρ_(vol-weighted)) for the fluid transfer by dividing the             accumulated volume-density product ({dot over (ν)}ρ_(accum))             by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, determining the non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer occurs after a fluid transfer end signal is received.

Preferably, the processing system is further configured to compare a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold and determine the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.

Preferably, the processing system is further configured to total aerated time portions to obtain an aerated transfer time (t_(aerated)), total all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(i)) and determine an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).

Preferably, the processing system is further configured to total aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)), total all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)), and determine an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).

Preferably, the processing system is further configured to obtain a temperature (T_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then add a volume-temperature product ({dot over (ν)}T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)), and determine a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the processing system is further configured to convert the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

Preferably, the processing system is further configured to obtain a fluid viscosity (μ_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then add a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)), and determine a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the processing system is further configured to receive the viscosity (μ_(i)) from an external source.

Preferably, the processing system is further configured to measure the viscosity with the vibratory flowmeter.

Preferably, the processing system is further configured to obtain a fluid viscosity (μ_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then add a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)), determine a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)), and convert the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

In one aspect of the invention, meter electronics for quantifying a fluid being transferred comprises:

-   -   an interface configured to communicate with a flowmeter assembly         of a vibratory flowmeter and receive a vibrational response; and     -   a processing system coupled to the interface and configured to:         -   measure a volume flow ({dot over (ν)}_(i)) for a             predetermined time portion (t_(i)) of the fluid transfer;         -   obtain a fluid viscosity (μ_(i)) for the predetermined time             portion (t_(i));         -   determine if the fluid transfer is non-aerated during the             predetermined time portion (t_(i));         -   if the predetermined time portion (t_(i)) is non-aerated,             then add a volume-viscosity product ({dot over             (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product             ({dot over (ν)}μ_(accum)) and add the volume flow ({dot over             (ν)}_(i)) to an accumulated volume flow ({dot over             (ν)}_(accum)); and         -   determine a non-aerated volume-weighted viscosity             (μ_(weighted)) for the fluid transfer by dividing the             accumulated volume-viscosity product ({dot over             (ν)}μ_(accum)) by the accumulated volume flow ({dot over             (ν)}_(accum)).

Preferably, determining the non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer occurs after a fluid transfer end signal is received.

Preferably, the processing system is further configured to compare a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold and determine the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.

Preferably, the processing system is further configured to total aerated time portions to obtain an aerated transfer time (t_(aerated)) total all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(total)), and determine an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).

Preferably, the processing system is further configured to total aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)), total all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)), and determine an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}^(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).

Preferably, the processing system is further configured to obtain a temperature (T_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then add a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)), and determine a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the processing system is further configured to convert the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

Preferably, the processing system is further configured to measure a density (ρ_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then add a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρaccum) and determine a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the processing system is further configured to measure a density (ρ_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then add a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)) 9 determine a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)), and convert the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

Preferably, the processing system is further configured to receive the viscosity (μ_(i)) from an external source.

Preferably, the processing system is further configured to measure the viscosity (μ_(i)) with the vibratory flowmeter.

In one aspect of the invention, a fluid quantification method for a fluid being transferred comprises:

-   -   measuring a volume flow ({dot over (ν)}_(i)) for a predetermined         time portion (t_(i)) of the fluid transfer, with the measuring         being performed by a vibratory flowmeter;     -   determining if the fluid transfer is non-aerated during the         predetermined time portion (t_(i));     -   if the predetermined time portion (t_(i)) is non-aerated, then         adding a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an         accumulated volume-density product ({dot over (ν)}ρ_(accum)) and         adding the volume flow ({dot over (ν)}_(i)) to an accumulated         volume flow ({dot over (ν)}_(accum)); and     -   determining a non-aerated volume-weighted density         (ρ_(vol-weighted)) for the fluid transfer by dividing the         accumulated volume-density product ({dot over (ν)}ρ_(accum)) by         the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, determining the non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer occurs after a fluid transfer end signal is received.

Preferably, determining if the fluid transfer is non-aerated comprises comparing a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold and determining the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.

Preferably, the method further comprises totaling aerated time portions to obtain an aerated transfer time (t_(aerated)), totaling all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(total)) and determining an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).

Preferably, the method further comprises totaling aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)), totaling all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)), and determining an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).

Preferably, the method further comprises obtaining a temperature (T_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)), and determining a non-aerated volume-weighted temperature (T_(vol-weighted) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T) _(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the method further comprises converting the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

Preferably, the method further comprises obtaining a viscosity (μ_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)), and determining a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, obtaining the viscosity (μ_(i)) comprises receiving the viscosity (μ_(i)) from an external source.

Preferably, obtaining the viscosity (μ_(i)) comprises measuring the viscosity (μ_(i)) with the vibratory flowmeter.

Preferably, the method further comprises obtaining a viscosity (μ_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)), determining a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)), and converting the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

In one aspect of the invention, a fluid quantification method for a fluid being transferred comprises:

-   -   measuring a volume flow ({dot over (ν)}_(i)) for a predetermined         time portion (t_(i)) of the fluid transfer, with the measuring         being performed by a vibratory flowmeter;     -   obtaining a viscosity (μ_(i)) for the predetermined time portion         (t_(i));     -   determining if the fluid transfer is non-aerated during the         predetermined time portion (t_(i));     -   if the predetermined time portion (t_(i)) is non-aerated, then         adding a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to         an accumulated volume-viscosity product ({dot over         (ν)}μ_(accum)) and adding the volume flow ({dot over (ν)}_(i))         to an accumulated volume flow ({dot over (ν)}_(accum)); and     -   determining a non-aerated volume-weighted viscosity         (μ_(vol-weighted)) for the fluid transfer by dividing the         accumulated volume-viscosity product ({dot over (ν)}μ_(accum))         by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, determining the non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer occurs after a fluid transfer end signal is received.

Preferably, determining if the fluid transfer is non-aerated comprises comparing a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold and determining the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.

Preferably, the method further comprises totaling aerated time portions to obtain an aerated transfer time (t_(aerated)), totaling all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(total)), and determining an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time a (t_(total)).

Preferably, the method further comprises totaling aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)) totaling all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)), and determining an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).

Preferably, the method further comprises obtaining a temperature (T_(i)) for the predetermined time portion (t_(i)), if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)), and determining a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the method further comprises converting the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

Preferably, the method further comprises measuring a density (ρ_(i)) for the predetermined time portion (t_(i)), with the measuring being performed by the vibratory flowmeter, if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)), and determining a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).

Preferably, the method further comprises measuring a density (ρ_(i)) for the predetermined time portion (t_(i)), with the measuring being performed by the vibratory flowmeter, if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)), determining a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)), and converting the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).

Preferably, obtaining the viscosity (μ_(i)) comprises receiving the viscosity (μ_(i)) from an external source.

Preferably, obtaining the viscosity (μ_(i)) comprises measuring the viscosity (μ_(i)) with the vibratory flowmeter.

DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a vibratory flowmeter according to the invention.

FIG. 2 shows a bunker fuel transfer according to the invention.

FIG. 3 is a block diagram of the meter electronics of the vibratory flowmeter according to the invention.

FIG. 4 is a flowchart of a fluid quantification method for a fluid being transferred according to the invention.

DETAILED DESCRIPTION OF THE INVENTION

FIGS. 1-4 and the following description depict specific examples to teach those skilled in the art how to make and use the best mode of the invention. For the purpose of teaching inventive principles, some conventional aspects have been simplified or omitted. Those skilled in the art will appreciate variations from these examples that fall within the scope of the invention. Those skilled in the art will appreciate that the features described below can be combined in various ways to form multiple variations of the invention. As a result, the invention is not limited to the specific examples described below, but only by the claims and their equivalents.

FIG. 1 shows a vibratory flowmeter 5 according to the invention. The vibratory flowmeter 5 comprises a flowmeter assembly 10 and meter electronics 20. The meter electronics 20 is connected to the meter assembly 10 via leads 100 and is configured to provide measurements of one or more of a density, mass flow rate, volume flow rate, totalized mass flow, temperature, or other measurements or information over a communication path 26. The vibratory flowmeter 5 can comprise a Coriolis mass flow meter. In addition, it should be recognized that the vibratory flowmeter 5 can alternatively comprise a vibratory densitometer. It should be apparent to those skilled in the art that the vibratory flowmeter can comprise any manner of vibratory flowmeter, regardless of the number of drivers, pick-off sensors, flow conduits, or the operating mode of vibration.

The flowmeter assembly 10 includes a pair of flanges 101 and 101′, manifolds 102 and 102′, a driver 104, pick-off sensors 105 and 105′, and flow conduits 103A and 103B. The driver 104 and the pick-off sensors 105 and 105′ are connected to the flow conduits 103A and 103B.

The flanges 101 and 101′ are affixed to the manifolds 102 and 102′. The manifolds 102 and 102′ can be affixed to opposite ends of a spacer 106 in some embodiments. The spacer 106 maintains the spacing between the manifolds 102 and 102′ in order to prevent pipeline forces from being transmitted to flow conduits 103A and 103B. When the flowmeter assembly 10 is inserted into a pipeline (not shown) which carries the flow fluid being measured, the flow fluid enters the flowmeter assembly 10 through the flange 101, passes through the inlet manifold 102 where the total amount of flow fluid is directed to enter the flow conduits 103A and 103B, flows through the flow conduits 103A and 103B and back into the outlet manifold 102′, where it exits the meter assembly 10 through the flange 101′.

The flow fluid can comprise a liquid. The flow fluid can comprise a gas. The flow fluid can comprise a multi-phase fluid, such as a liquid including entrained gases and/or entrained solids.

The flow conduits 103A and 103B are selected and appropriately mounted to the inlet manifold 102 and to the outlet manifold 102′ so as to have substantially the same mass distribution, moments of inertia, and elastic modules about the bending axes W-W and W′-W′ respectively. The flow conduits 103A and 103B extend outwardly from the manifolds 102 and 102′ in an essentially parallel fashion.

The flow conduits 103A and 103B are driven by the driver 104 in opposite directions about the respective bending axes W and W′ and at what is termed the first out of phase bending mode of the vibratory flowmeter 5. The driver 104 may comprise one of many well known arrangements, such as a magnet mounted to the flow conduit 103A and an opposing coil mounted to flow conduit 103B. An alternating current is passed through the opposing coil to cause both conduits to oscillate. A suitable drive signal is applied by the meter electronics 20 to the driver 104 via the lead 110. Other driver devices are contemplated and are within the scope of the description and claims.

The meter electronics 20 receives sensor signals on the leads 111 and 111′, respectively. The meter electronics 20 produces a drive signal on the lead 110 which causes the driver 104 to oscillate the flow conduits 103A and 103B. Other sensor devices are contemplated and are within the scope of the description and claims.

The meter electronics 20 processes the left and right velocity signals from the pick-off sensors 105 and 105′ in order to compute a flow rate, among other things. The communication path 26 provides an input and an output means that allows the meter electronics 20 to interface with an operator or with other electronic systems. The description of FIG. 1 is provided merely as an example of the operation of a Coriolis flowmeter and is not intended to limit the teaching of the present invention.

The meter electronics 20 in one embodiment is configured to vibrate the flowtubes 103A and 103B. The vibration is performed by the driver 104. The meter electronics 20 further receives resulting vibrational signals from the pickoff sensors 105 and 105′. The vibrational signals comprise a vibrational response of the flowtubes 103A and 103B. The meter electronics 20 processes the vibrational response and determines a response frequency and/or phase difference. The meter electronics 20 processes the vibrational response and determines one or more flow measurements, including a mass flow rate and/or density of the flow fluid. Other vibrational response characteristics and/or flow measurements are contemplated and are within the scope of the description and claims.

In one embodiment, the flowtubes 103A and 103B comprise substantially U-shaped flowtubes, as shown. Alternatively, in other embodiments, the flowtubes can comprise substantially straight flowtubes. Additional flowmeter shapes and/or configurations can be used and are within the scope of the description and claims.

FIG. 2 shows a bunker fuel transfer according to the invention. The bunker fuel transfer comprises a transfer and metering of fuel from a first tank 3 to a second tank 4. In the figure, the bunker fuel is being transferred from a fueling barge 1 to a ship 2, but it should be understood that the transfer can comprise any fuel transfer. It should be understood that the transfer can comprise any fluid transfer. In addition, it should be understood that a flow meter may be installed on the barge, the ship, or in both locations.

Bunkering is generally known as the practice of storing and transferring marine fuel oils, which have come to be known as bunker fuels. For ship fueling, large amounts of fuel may be temporarily stored in a barge or other container for the purpose of transferring fuel from shore to a ship. A bunker may be located on a dock or other port facility, or may be carried by a barge or other refueling vehicle. During bunkering, the fuel measurement usually comprises an empty-full-empty batching process, allowing gas to become entrained in the fuel.

Bunker fuel comprises a relatively heavy petroleum derivative that is used in heating or in large industrial and/or marine engines. There are multiple grades of fuel that may comprise a bunker fuel. Bunker fuel is generally heavier and more viscous than gasoline or diesel.

Marine fuel costs represent a major portion of a ship's operating cost. With increasing oil prices and increasing conservation efforts, careful fuel management has become vital for environmental and financial reasons.

Blending a combination of fuel components in exact proportions for use in the marine industry typically involves the combination of heavy fuel oil (HFO or Bunker C) and an intermediate fuel oil (IFO), which is available in a range of viscosities and sulfur contents. The International Standard Organization (ISO) has specifications for marine fuels for worldwide use onboard ships. Precision blending is very important due to increases in fuel prices, the danger of engine wear from different fuel grade burn temperatures, and the limitations associated with shipboard storage. Preload blending guarantees that a vessel receives fuel with optimal properties for the intended use of specific onboard engines, leading to reductions in engine NOx and Sox emissions. The blending of fuels pre-bunkering has also been increasing in popularity in the marine industry due to regulations discouraging the mixing of fuels onboard the vessel. Blending operations often require high accuracy and a meter which is relatively immune to entrained gas, making mass flow meters ideal.

FIG. 3 is a block diagram of the meter electronics 20 of the vibratory flowmeter 5 according to the invention. In operation, the vibratory flowmeter 5 is used to quantify a fluid during a fluid transfer. The fluid can include a fuel. The vibratory flowmeter 5 can be used to measure a volume flow ({dot over (ν)}_(i)) value and/or volume flow total ({dot over (ν)}_(tot)) of a fluid transfer. In some embodiments the fluid transfer is measured over a series of predetermined time portions (t_(i)). The predetermined time portions (t_(i)) can be uniform or non-uniform in duration. The predetermined time portions (t_(i)) can be chosen to generate a plurality of measurements during the transfer. The length of a predetermined time portion (t_(i)) can be chosen so as to capture accurate and representative values.

The vibratory flowmeter 5 generates a vibrational response. The vibrational response is received and processed by the meter electronics 20 to generate one or more fluid quantification values. The values can be monitored, recorded, and totaled.

The meter electronics 20 includes an interface 201, a processing system 203 in communication with the interface 201, and a storage system 204 in communication with the processing system 203. Although these components are shown as distinct blocks, it should be understood that the meter electronics 20 can be comprised of various combinations of integrated and/or discrete components.

The interface 201 is configured to communicate with the flowmeter assembly 10 of the vibratory flowmeter 5. The interface 201 may be configured to couple to the leads 100 (see FIG. 1) and exchange signals with the driver 104 and the pickoff sensors 105 and 105′. The interface 201 may be further configured to communicate over the communication path 26, such as to external devices.

The processing system 203 can comprise any manner of processing system. The processing system 203 is configured to retrieve and execute stored routines 205 in order to operate the vibratory flowmeter 5. The storage system 204 can store routines including a flowmeter routine 205, a volume-weighted density routine 209, a volume-weighted viscosity routine 210, a volume-weighted temperature routine 211, and an aeration detect routine 213. Other measurement/processing routines are contemplated and are within the scope of the description and claims. The storage system 204 can store measurements, received values, working values, and other information. In some embodiments, the storage system stores a volume flow ({dot over (ν)}_(i)) 221, a density (ρ) 222, a viscosity (μ) 223, a temperature (T) 224, a volume-density product ({dot over (ν)}_(i)ρ_(i)) 234, a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) 235, a volume-temperature product ({dot over (ν)}_(i)T_(i)) 236, a volume-weighted density (ρ_(vol-weighted)) 241, a volume-weighted viscosity (μ_(vol-weighted)) 242, a volume-weighted temperature (T_(vol-weighted)) 243, an aeration threshold 244, and an aeration fraction 248.

The flowmeter routine 205 can produce and store fluid quantifications and flow measurements. These values can comprise substantially instantaneous measurement values or can comprise totaled or accumulated values. For example, the flowmeter routine 205 can generate volume flow measurements and store them in the volume flow ({dot over (ν)}_(i)) storage 221. The flowmeter routine 205 can generate density measurements and store them in the density (ρ) storage 222. The volume flow ({dot over (ν)}_(i)) and density (ρ_(i)) values are determined from the vibrational response, as previously discussed and as known in the art. The volume flow ({dot over (ν)}_(i)) can comprise a substantially instantaneous volume flow rate value, can comprise a volume flow rate sample, can comprise an averaged volume flow rate over the time portion (t_(i)), or can comprise an accumulated volume flow rate over the time portion (t_(i)). In addition, other volume flow quantifications are contemplated and are within the scope of the description and claims.

The flowmeter routine 205 can generate or receive temperature measurements and store them in the temperature (T) storage 224. The flowmeter routine 205 can generate or receive viscosity measurements and store them in the viscosity (μ) storage 223.

The vibratory flowmeter 5 in some embodiments can be configured to perform viscosity measurements. One such vibratory flowmeter is disclosed in co-pending U.S. Patent Publication No. 2008/0184813, assigned to the current assignee. The 2008/0184813 patent publication is herein incorporated by reference.

Alternatively, in other embodiments the flowmeter routine 205 receives viscosity values from an external source and stores the received values in the viscosity (μ) storage 223. The external source can be a separate viscosity meter, from an operator input, from a stored value, or other sources.

In some embodiments, the meter electronics 20 is configured to execute the volume-weighted density routine 209. The volume-weighted density routine 209 is configured to measure a volume flow ({dot over (ν)}_(i)) and a density (ρ_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with the volume flow ({dot over (ν)}_(i)) and density (ρ_(i)) measurements being performed by the vibratory flowmeter 5 through vibration of the flowmeter assembly 10, determine a volume-density product ({dot over (ν)}_(i)ρ_(i)) for the predetermined time portion (t_(i)), accumulate the volume flow ({dot over (ν)}_(i)) and the volume-density product ({dot over (ν)}_(i)ρ_(i)) values into an accumulated volume flow ({dot over (ν)}_(accum)) and an accumulated volume-density product ({dot over (m)}ρ_(accum)) for all non-aerated time portions of the fluid transfer, and determine a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)). The resulting non-aerated volume-weighted density (ρ_(vol-weighted)) can be stored in the volume-weighted density storage 238. This processing may additionally generate the accumulated volume-density product ({dot over (ν)}ρ_(accum)) that can be stored in the volume-density product storage 234. The processing therefore generates the volume-weighted density (ρ_(vol-weighted)) according to:

$\begin{matrix} {\rho_{{vol} - {weighted}} = {\frac{\sum{{\overset{.}{v}}_{i}\rho_{i}}}{\sum{\overset{.}{v}}_{i}} = \frac{\left( {{\overset{.}{v}}_{1}\rho_{1}} \right) + \left( {{\overset{.}{v}}_{2}\rho_{2}} \right) + \left( {{\overset{.}{v}}_{3}\rho_{3}} \right) + {\ldots \mspace{14mu} \left( {{\overset{.}{v}}_{N}\rho_{N}} \right)}}{{\overset{.}{v}}_{1} + {\overset{.}{v}}_{2} + {\overset{.}{v}}_{3} + {\ldots \mspace{14mu} {\overset{.}{v}}_{N}}}}} & (1) \end{matrix}$

The volume-weighted density (ρ_(vol-weighted)) can be similar to an average density, representing the density of the fluid transfer, but where only the density measurements for non-aerated time portions are taken into account. Density measurements from aerated time portions can be left out, as such aerated density measurements are likely to be inaccurate and unreliable. Also, the volume-weighted quantity is superior to a basic time-averaged or periodically sampled density because it represents a true average quantity with respect to the volume of fuel purchased. A time-averaged algorithm would be sensitive to variations in flow rate, and a periodic sampling system is sensitive to variation in density due to stratification of tanks.

When the vibratory flowmeter 5 is operated as a densitometer or viscosity meter, aeration in the fluid can affect and degrade density and viscosity measurements (and possibly temperature measurements). For this reason, the volume-density, volume-viscosity, and/or volume-temperature products are not totaled during aeration time periods.

In some embodiments, the meter electronics 20 is configured to execute the volume-weighted viscosity routine 210. The volume-weighted viscosity routine 210 is configured to measure a volume flow ({dot over (ν)}_(i)) and obtain a fluid viscosity (μ_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with at least the volume flow ({dot over (ν)}_(i)) measurement being performed by the vibratory flowmeter 5 through vibration of the flowmeter assembly 10, determine a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) for the predetermined time portion (t_(i)), accumulated the volume flow ({dot over (ν)}_(i)) and the volume-viscosity product ({dot over (ν)}_(i)μ_(i)) values into an accumulated volume flow ({dot over (ν)}_(accum)) and an accumulated volume-viscosity product ({dot over (ν)}_(accum)) for all non-aerated predetermined time portions of the fluid transfer, and determine a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)). This non-aerated volume-weighted viscosity (μ_(vol-weighted)) can be stored in the volume-weighted viscosity storage 239. The viscosity measurements can reflect quality or grade of the flow fluid as it is transferred. This processing may additionally generate the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) that can be stored in the volume-viscosity product storage 235. The processing therefore generates the volume-weighted viscosity (μ_(vol-weighted)) according to:

$\begin{matrix} {\mu_{{vol} - {weighted}} = \frac{\sum{{\overset{.}{v}}_{i}\mu_{i}}}{\sum{\overset{.}{v}}_{i}}} & (2) \end{matrix}$

The volume-weighted viscosity (μ_(vol-weighted)) can be similar to an average viscosity, representing the overall viscosity of the fluid transfer, but where only the viscosity measurements for non-aerated time portions are taken into account. Viscosity measurements from aerated time portions can be left out, as such aerated viscosity measurements are likely to be inaccurate and unreliable.

In some embodiments, the meter electronics 20 is configured to execute the volume-weighted temperature routine 211. The volume-weighted temperature routine 211 is configured to measure a volume flow ({dot over (ν)}_(i)) and a temperature (T_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with at least the volume flow ({dot over (ν)}_(i)) measurement being performed by the vibratory flowmeter 5 through vibration of the flowmeter assembly 10, determine a volume-temperature product ({dot over (ν)}_(i)T_(i)) for the predetermined time portion (t_(i)), accumulate the volume flow ({dot over (ν)}_(i)) and the volume-temperature product ({dot over (ν)}_(i)T_(i)) values into an accumulated volume flow ({dot over (ν)}_(accum)) and an accumulated volume-temperature product ({dot over (ν)}T_(accum)) for all non-aerated predetermined time portions of the fluid transfer, and determine a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)). This non-aerated volume-weighted temperature (T_(vol-weighted)) can be stored in the volume-weighted temperature storage 243. This processing may additionally generate the accumulated volume-temperature product ({dot over (ν)}T_(accum)) that can be stored in the volume-temperature product storage 236. The processing therefore generates the volume-weighted temperature (T_(vol-weighted)) according to:

$\begin{matrix} {T_{{vol} - {weighted}} = \frac{\sum{{\overset{.}{v}}_{i}T_{i}}}{\sum{\overset{.}{v}}_{i}}} & (3) \end{matrix}$

The volume-weighted temperature (T_(vol-weighted)) can be similar to an average temperature, representing the overall temperature of the fluid transfer, but where only the temperature measurements for non-aerated time portions are taken into account. Temperature measurements from aerated time portions can be left out.

In some embodiments, the meter electronics 20 is configured to execute the aeration detect routine 213. The aeration detect routine 213 processes flowmeter values in order to detect substantial aeration in a fluid transfer. The detection can be substantially continuously performed during the fluid transfer, such as during each time portion (t_(i)). The aeration detect routine 213 is configured to vibrate the flowmeter assembly 10 of a vibratory flowmeter 5 during the fluid transfer, compare a vibrational response of each time portion (t_(i)) of the fluid transfer to an aeration threshold 244, and determine that a time portion (t_(i)) is an aerated time portion where the vibrational response does not exceed the aeration threshold 244. The aeration threshold 244 comprises a range above which the fluid is considered to be non-aerated or is considered to be minimally (e.g., acceptably) aerated.

In some embodiments, the aeration threshold 244 comprises a drive gain threshold. A current drive gain is compared to the drive gain threshold. The drive gain quantifies the amount of response per a given drive input. The drive gain can comprise a pickoff response amplitude divided by a driver vibration amplitude, where the values can be instantaneous or can be averaged or integrated over the predetermined time portion (t_(i)). The drive gain will typically be stable at around ten to thirty percent for a pure liquid and will increase dramatically as entrained gas enters the flow fluid. If the drive gain fails the aeration threshold 244, then the meter electronics 20 can determine that the fluid is aerated. For example, the drive gain may fail the aeration threshold 244 if the drive gain exceeds the aeration threshold 244, as the drive gain may increase in value with the presence of aeration.

In some embodiments, the aeration threshold 244 comprises a minimum response amplitude. The minimum response amplitude relates to the vibrational response amplitude generated by a pickoff sensor 105 or 105′. It is known that the pickoff amplitude will decrease with entrained air in the flow fluid. If the pickoff amplitude fails the aeration threshold 244, then the meter electronics 20 can determine that the fluid is aerated. For example, the pickoff amplitude may fail the aeration threshold 244 if the pickoff amplitude is less than the aeration threshold 244, as the pickoff amplitude may decrease in value with the presence of aeration.

In some embodiments, the aeration threshold 244 comprises a minimum density threshold for the fluid. The predetermined aerated density threshold can comprise a level of acceptable aeration in the fluid, such as for air in a bunker fuel, for example. Consequently, the measured density (ρ_(i)) can be compared to the aeration threshold 244. If the measured density (ρ_(i)) fails the aeration threshold 244, then the meter electronics 20 can determine that the fluid is aerated. For example, the measured density (ρ_(i)) may fail the aeration threshold 244 if the measured density (ρ_(i)) is less than the aeration threshold 244, as the fluid density will decrease with aeration.

Aeration may be detected based on a damping estimate, found from the shape of a frequency response function for the drive mode of the flow meter. For example, the width of the drive mode peak in a frequency response graph can help determine whether gas is present. A wider peak indicates higher damping and therefore the presence of gas. Consequently, the damping estimate/drive mode peak width can be compared to the aeration threshold 244. If the damping estimate/drive mode peak width fails the aeration threshold 244, then the meter electronics 20 can determine that the fluid is aerated. For example, the damping estimate/drive mode peak width may fail the aeration threshold 244 if the damping estimate/drive mode peak width exceeds the aeration threshold 244, as the damping may increase in the presence of aeration.

Alternatively, a 3-dB down point may be used, which refers to the width of the peak at a specific off-resonant frequency which corresponds to a specific vibrational response amplitude. Consequently, the 3-dB down point can be compared to the aeration threshold 244. If the 3-dB down point fails the aeration threshold 244, then the meter electronics 20 can determine that the fluid is aerated. For example, the 3-dB down point may fail the aeration threshold 244 if the 3-dB down point is less than the aeration threshold 244, as the magnitude of the 3-dB down point may drop in the presence of aeration.

Flow noise increases with increasing aeration, and therefore it is also possible to detect the presence of gas by analysis of the standard deviations in the mass flow, density, or other measurements from a Coriolis meter. A high standard deviation can indicate an increase in aeration. Consequently, the standard deviation in the flow measurements due to flow noise (or other disturbances), can be compared to the aeration threshold 244. If the standard deviation fails the aeration threshold 244, then the meter electronics 20 can determine that the fluid is aerated. For example, the standard deviation may fail the aeration threshold 244 if the standard deviation exceeds the aeration threshold 244.

An aeration fraction is determined in some embodiments and stored in the aeration fraction storage 248. The aeration fraction denotes an amount of the fluid transfer that is aerated, and may serve as another fluid quality indicator. In some embodiments, the volume flow ({dot over (ν)}_(i)) values can be totaled for time portions (t_(i)) of detected aeration. The aerated volume ({dot over (ν)}_(aerated)) can be divided by the total volume ({dot over (ν)}_(tot)) to provide an aeration volume fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)), i.e., a measure of the amount of transferred volume that is aerated. Alternatively, the aerated time portions can be totaled and then divided by a total transfer time to generate an aeration time fraction, for example. Other aeration fraction quantifications are contemplated and are within the scope of the description and claims.

The volume-weighted and/or non-aerated measurements will provide better and more accurate measurements of a fluid transfer than strictly volumetric flow measurements or non-aeration discriminating mass flow measurements. The volume-weighted measurements and the accumulated non-aerated accumulated volume flow ({dot over (ν)}_(accum)) will provide a better and more accurate measurement of the energy content of the transferred fuel, in contrast to a non-aeration discriminating volume measurement of the prior art. A prior art non-aeration discriminating volume measurement will not take into account any aeration of the fuel.

FIG. 4 is a flowchart 400 of a fluid quantification method for a fluid being transferred according to the invention. The method comprises taking measurements at predetermined time portions (t_(i)). In step 401, volume flow ({dot over (ν)}_(i)), density (ρ_(i)), and temperature (T_(i)) values are measured for a current predetermined time portion (t_(i)). The volume flow ({dot over (ν)}_(i)) and density (ρ_(i)) values are determined from the vibrational response, as previously discussed and as known in the art. For example, a volume flow rate ({dot over (ν)}) can comprise

$\overset{.}{v} = {\frac{\overset{.}{m}}{\rho}.}$

Step 401 can comprise a start of the fluid transfer process. In addition, step 401 can comprise an iterative measurement step obtained at any time during the fluid transfer process.

The predetermined time portions (t_(i)) can be of any desired length that will adequately characterize the fluid being transferred. The predetermined time portions (t_(i)) are preferably substantially uniform in length, but do not have to be uniform or consistent.

In step 402, a viscosity (μ_(i)) value is obtained. The viscosity (μ_(i)) may be measured. Alternatively, the viscosity (μ_(i)) may be received from an external measurement source such as a viscometer or a differential pressure measurement across the flow meter or a length of pipe.

In step 403, a volume-density product ({dot over (ν)}_(i)ρ_(i)) is generated from the volume flow ({dot over (ν)}_(i)) and density ({dot over (ν)}_(i)) measurements. A volume-viscosity product ({dot over (ν)}_(i)μ_(i)) is generated from the volume flow ({dot over (ν)}_(i)) and viscosity (μ_(i)) measurements. A volume-temperature product ({dot over (ν)}_(i)T_(i)) is generated from the volume flow ({dot over (ν)}_(i)) and temperature (T_(i)) measurements. The volume flow ({dot over (ν)}_(i)), density (ρ_(i)), volume-density product ({dot over (ν)}_(i)ρ_(i)), volume-viscosity product ({dot over (ν)}_(i)μ_(i)), and volume-temperature product ({dot over (ν)}_(i)T_(i)) values can be added into previous values and therefore accumulated. Alternatively, the volume flow ({dot over (ν)}_(i)), density (ρ_(i)), volume-density product ({dot over (ν)}_(i)ρ_(i)), volume-viscosity product ({dot over (ν)}_(i)μ_(i)), and volume-temperature product ({dot over (ν)}_(i)T_(i)) values can be selectively accumulated, as discussed below.

In step 404, if the fluid transfer during the current time portion (t_(i)) is not substantially aerated, then the method proceeds to step 405. Otherwise, where the fluid is determined to be aerated, the method branches back to step 401 and the volume flow ({dot over (ν)}_(i)), density (μ_(i)), volume-density product ({dot over (ν)}_(i)ρ_(i)), volume-viscosity product ({dot over (ν)}_(i)μ_(i)), and volume-temperature product ({dot over (ν)}_(i)T_(i)) values are not accumulated. Instead, new values are obtained upon a next predetermined time portion (t_(i+1)). However, it should be understood that the measurements may be recorded in some manner and therefore may be available for other use, if desired.

In an alternative embodiment, if the fluid is determined to be substantially aerated, then the measurements may simply not be performed during the current time portion (t_(i)). In this embodiment, step 403 could be performed as a first step and would loop back on itself if aeration were detected, with the measuring/obtaining of steps 401, 402, and 403 being skipped in the event of aeration detection. The volume flow ({dot over (ν)}_(i)), density (ρ_(i)), volume-density product ({dot over (ν)}_(i)ρ_(i)), volume-viscosity product ({dot over (ν)}_(i)μ_(i)), and volume-temperature product ({dot over (ν)}_(i)T_(i)) values could simply be neglected, saving processing time.

In step 405, the volume flow ({dot over (ν)}_(i)), density (ρ_(i)), volume-density product ({dot over (ν)}_(i)ρ_(i)), volume-viscosity product ({dot over (ν)}_(i)μ_(i)), and volume-temperature product ({dot over (ν)}_(i)T_(i)) values are accumulated, added into respective running totals for the fluid transfer. The accumulated volume flow ({dot over (ν)}_(accum)), when the fluid transfer is completed, in some embodiments may represent the volume of non-aerated fluid transferred. The total density (ρ_(tot)) can be processed when the fluid transfer is completed in order to determine a mean, average, weighted, or volume-weighted density of fluid transferred. The volume-weighted density is desirable because the fluid being transferred may not be uniform in nature. A measure or quantification of a density for the fluid transfer as a whole may be more useful than an averaged density. Moreover, elimination of the volume values and/or density values of an aerated portion or portions of the fluid transfer may avoid skewing or mis-characterizing the fluid. If the fluid is highly aerated, the density measurements will be affected. Not only will the quantity of delivered fluid decrease when aerated, but also the vibratory flowmeter will generate inaccurate density measurements due to the aeration.

In step 406, if the method is done, i.e., the fluid transfer is complete, then the method proceeds to step 407. Otherwise, where the fluid transfer is not complete, the method branches back to step 401. In this manner, the measurement process is iteratively performed over the entire fluid transfer process.

The fluid transfer may be complete when a fluid transfer end signal is received in the vibratory flowmeter. The fluid transfer end signal may be received from an operator or from another device. Alternatively, the vibratory flowmeter can autonomously determine an end of the fluid transfer and generate the fluid transfer end signal. For example, the vibratory flowmeter can autonomously generate the fluid transfer end signal in some embodiments if the volume flow drops below a predetermined transfer threshold for more than a predetermined time period.

In step 407, where the fluid transfer is complete, values are determined for the fluid transfer as a whole. A non-aerated volume-weighted density (ρ_(vol-weighted)) can be determined for the fluid transfer. The non-aerated volume-weighted density (ρ_(vol-weighted)) comprises the accumulated volume-density product ({dot over (ν)}ρ_(accum)) divided by the accumulated volume flow ({dot over (ν)}_(accum)). This can be represented by equation (1), above. The non-aerated volume-weighted density (ρ_(vol-weighted)) in some embodiments comprises a density generated from only the non-aerated fluid portions. The non-aerated volume-weighted density (ρ_(vol-weighted)) therefore comprises a more accurate density for the fluid when there is any aeration in the fluid. If the fluid is not determined to be aerated, then the volume-weighted density (ρ_(vol-weighted)) will substantially match an average density for the entire fluid transfer.

A non-aerated volume-weighted viscosity (μ_(vol-weighted)) can be determined for the fluid transfer. The non-aerated volume-weighted viscosity (μ_(vol-weighted)) comprises the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) divided by the accumulated volume flow ({dot over (ν)}_(accum)). This can be represented by equation (2), above. The non-aerated volume-weighted viscosity (μ_(vol-weighted)) in some embodiments comprises a viscosity generated from only the non-aerated fluid portions.

A non-aerated volume-weighted temperature (T_(vol-weighted)) can be determined for the fluid transfer. The non-aerated volume-weighted temperature (T_(vol-weighted)) comprises the accumulated volume-temperature product ({dot over (ν)}T_(accum)) divided by the accumulated volume flow ({dot over (ν)}_(accum)). This can be represented by equation (3), above. The non-aerated volume-weighted temperature (T_(vol-weighted)) in some embodiments comprises a temperature generated from only the non-aerated fluid portions.

In step 408, the non-aerated volume-weighted temperature (T_(vol-weighted)) can be used to perform temperature compensation. For example, the computed volume-weighted density (ρ_(vol-weighted)) can be temperature-compensated to generate a volume-weighted density value for a standard temperature. In addition or alternatively, the volume-weighted viscosity (μ_(vol-weighted)) can also be temperature compensated. This can be done using known methods, such as by using a table or tables available from the American Petroleum Institute (API) for a bunker fuel or other fuel product, wherein the volume-weighted density and the volume-weighted temperature are used to produce a corresponding standard density value. The volume-weighted viscosity and the volume-weighted temperature are used to produce a corresponding standard viscosity value. Such a table can be used to translate a value at a given temperature to the value at a reference temperature. The reference temperature value is useful for comparing the computed volume-weighted density and/or the computed volume-weighted viscosity to a standard value in order to judge the quality of the bunker fuel or other fluid. In this manner, the relative quality of the bunker fuel can be assessed, where if the computed volume-weighted density/viscosity deviates significantly from a standard value, then the quality of the bunker fuel (or other fluid) may be poor or unacceptable.

In step 409, an aeration fraction is determined, as previously discussed.

The meter electronics and fluid quantification method according to the invention can be employed according to any of the embodiments in order to provide several advantages, if desired. The meter electronics and fluid quantification method can generate a more reliable measurement of potentially aerated fluids. The meter electronics and fluid quantification method can generate a more accurate measurement of fluid properties, even if those properties vary throughout the batch. The resulting mass flow rate measurement is not affected by aeration or by ambient temperature or pressure. The resulting volume flow rate measurement is not affected by aeration. The meter electronics and fluid quantification method can generate a more accurate measurement of fluid transfer, wherein the measurement is not affected by aeration. The meter electronics and fluid quantification method can detect aeration in the fluid. The meter electronics and fluid quantification method can generate a better measure of energy content of a fuel. The meter electronics and fluid quantification method can cross-check a measured mass of the fluid transfer against a measured and/or estimated volume. The meter electronics and fluid quantification method can generate measurements and records of fluid variation during transfer. 

1. Meter electronics (20) for quantifying a fluid being transferred, the meter electronics (20) comprising an interface (201) configured to communicate with a flowmeter assembly of a vibratory flowmeter and receive a vibrational response and a processing system (203) coupled to the interface (201) and configured to measure a volume flow ({dot over (ν)}_(i)) and a density (ρ_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with the processing system (203) being characterized by being configured to: determine if the fluid transfer is non-aerated during the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)) and add the volume flow ({dot over (ν)}_(i)) to an accumulated volume flow ({dot over (ν)}_(accum)); and determine a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 2. The meter electronics (20) of claim 1, with determining the non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer occurring after a fluid transfer end signal is received.
 3. The meter electronics (20) of claim 1, with the processing system (203) being further configured to: compare a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold; and determine the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.
 4. The meter electronics (20) of claim 3, with the processing system (203) being further configured to: total aerated time portions to obtain an aerated transfer time (t_(aerated)); total all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(total)); and determine an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).
 5. The meter electronics (20) of claim 3, with the processing system (203) being further configured to: total aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)); total all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)); and determine an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).
 6. The meter electronics (20) of claim 1, with the processing system (203) being further configured to: obtain a temperature (T_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)); and determine a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 7. The meter electronics (20) of claim 6, with the processing system (203) being further configured to convert the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 8. The meter electronics (20) of claim 1, with the processing system (203) being further configured to: obtain a fluid viscosity (μ_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)); and determine a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 9. The meter electronics (20) of claim 8, with the processing system (203) being further configured to receive the viscosity (μ_(i)) from an external source.
 10. The meter electronics (20) of claim 8, with the processing system (203) being further configured to measure the viscosity (μ_(i)) with the vibratory flowmeter (5).
 11. The meter electronics (20) of claim 6, with the processing system (203) being further configured to: obtain a fluid viscosity (μ_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)); determine a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)); and convert the non-aerated volume-weighted viscosity ({dot over (ν)}_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 12. Meter electronics (20) for quantifying a fluid being transferred, the meter electronics (20) comprising an interface (201) configured to communicate with a flowmeter assembly of a vibratory flowmeter and receive a vibrational response and a processing system (203) coupled to the interface (201) and configured to measure a volume flow ({dot over (ν)}_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with the processing system (203) being characterized by being configured to: obtain a fluid viscosity (μ_(i)) for the predetermined time portion (t_(i)); determine if the fluid transfer is non-aerated during the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) and add the volume flow ({dot over (ν)}_(i)) to an accumulated volume flow ({dot over (ν)}_(accum)); and determine a non-aerated volume-weighted viscosity ({dot over (ν)}_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 13. The meter electronics (20) of claim 12, with determining the non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer occurring after a fluid transfer end signal is received.
 14. The meter electronics (20) of claim 12, with the processing system (203) being further configured to: compare a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold; and determine the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.
 15. The meter electronics (20) of claim 14, with the processing system (203) being further configured to: total aerated time portions to obtain an aerated transfer time (t_(aerated)); total all time portions (t_(i)) of the fluid transfer to obtain a total transfer time a (t_(total)); and determine an aeration fraction a (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).
 16. The meter electronics (20) of claim 14, with the processing system (203) being further configured to: total aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)); total all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)); and determine an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).
 17. The meter electronics (20) of claim 12, with the processing system (203) being further configured to: obtain a temperature (T_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)); and determine a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 18. The meter electronics (20) of claim 17, with the processing system (203) being further configured to convert the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 19. The meter electronics (20) of claim 12, with the processing system (203) being further configured to: measure a density (ρ_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)); and determine a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 20. The meter electronics (20) of claim 17, with the processing system (203) being further configured to: measure a density (ρ_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then add a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}_(accum)); determine a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)); and convert the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 21. The meter electronics (20) of claim 12, with the processing system (203) being further configured to receive the viscosity (μ_(i)) from an external source.
 22. The meter electronics (20) of claim 12, with the processing system (203) being further configured to measure the viscosity (μ_(i)) with the vibratory flowmeter (5).
 23. A fluid quantification method for a fluid being transferred, the method comprising measuring a volume flow ({dot over (ν)}_(i)) and a density (ρ_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with the measuring being performed by a vibratory flowmeter, with the method being characterized by: determining if the fluid transfer is non-aerated during the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)) and adding the volume flow ({dot over (ν)}_(i)) to an accumulated volume flow ({dot over (ν)}_(accum)); and determining a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 24. The method of claim 23, with determining the non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer occurring after a fluid transfer end signal is received.
 25. The method of claim 23, with determining if the fluid transfer is non-aerated comprising: comparing a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold; and determining the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.
 26. The method of claim 25, further comprising: totaling aerated time portions to obtain an aerated transfer time (t_(aerated)); totaling all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(total)); and determining an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).
 27. The method of claim 25, further comprising: totaling aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)); totaling all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)); and determine an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).
 28. The method of claim 23, further comprising: obtaining a temperature (T_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)); and determining a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 29. The method of claim 28, further comprising converting the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 30. The method of claim 23, further comprising: obtaining a viscosity (μ_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)); and determining a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 31. The method of claim 30, with obtaining the viscosity (μ_(i)) comprising receiving the viscosity (μ_(i)) from an external source.
 32. The method of claim 30, with obtaining the viscosity (μ_(i)) comprising measuring the viscosity (μ_(i)) with the vibratory flowmeter.
 33. The method of claim 28, further comprising: obtaining a viscosity (μ_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)); determining a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)); and converting the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 34. A fluid quantification method for a fluid being transferred, the method comprising measuring a volume flow ({dot over (ν)}_(i)) for a predetermined time portion (t_(i)) of the fluid transfer, with the measuring being performed by a vibratory flowmeter, with the method being characterized by: obtaining a viscosity (μ_(i)) for the predetermined time portion (t_(i)); determining if the fluid transfer is non-aerated during the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-viscosity product ({dot over (ν)}_(i)μ_(i)) to an accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) and adding the volume flow ({dot over (ν)}_(i)) to an accumulated volume flow ({dot over (ν)}_(accum)); and determining a non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-viscosity product ({dot over (ν)}μ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 35. The method of claim 34, with determining the non-aerated volume-weighted viscosity (μ_(vol-weighted)) for the fluid transfer occurring after a fluid transfer end signal is received.
 36. The method of claim 34, with determining if the fluid transfer is non-aerated comprising: comparing a vibrational response of the predetermined time portion (t_(i)) to a predetermined aeration threshold; and determining the predetermined time portion (t_(i)) to be aerated if the vibrational response fails the predetermined aeration threshold.
 37. The method of claim 36, further comprising: totaling aerated time portions to obtain an aerated transfer time (t_(aerated)); totaling all time portions (t_(i)) of the fluid transfer to obtain a total transfer time (t_(total)); and determining an aeration fraction (t_(aerated)/t_(total)) as the aerated transfer time (t_(aerated)) divided by the total transfer time (t_(total)).
 38. The method of claim 36, further comprising: totaling aerated volume flow ({dot over (ν)}_(i)) values to obtain an aerated volume flow ({dot over (ν)}_(aerated)); totaling all volume flow ({dot over (ν)}_(i)) values to obtain a total volume flow ({dot over (ν)}_(tot)); and determine an aeration fraction ({dot over (ν)}_(aerated)/{dot over (ν)}_(tot)) as the aerated volume flow ({dot over (ν)}_(aerated)) divided by the total volume flow ({dot over (ν)}_(tot)).
 39. The method of claim 34, further comprising: obtaining a temperature (T_(i)) for the predetermined time portion (t_(i)); if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-temperature product ({dot over (ν)}_(i)T_(i)) to an accumulated volume-temperature product ({dot over (ν)}T_(accum)); and determining a non-aerated volume-weighted temperature (T_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-temperature product ({dot over (ν)}T_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 40. The method of claim 39, further comprising converting the non-aerated volume-weighted viscosity (μ_(vol-weighted)) to a standard viscosity value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 41. The method of claim 34, further comprising: measuring a density (ρ_(i)) for the predetermined time portion (t_(i)), with the measuring being performed by the vibratory flowmeter; if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}_(accum)); and determining a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)).
 42. The method of claim 39, further comprising: measuring a density (ρ_(i)) for the predetermined time portion (t_(i)), with the measuring being performed by the vibratory flowmeter; if the predetermined time portion (t_(i)) is non-aerated, then adding a volume-density product ({dot over (ν)}_(i)ρ_(i)) to an accumulated volume-density product ({dot over (ν)}ρ_(accum)); determining a non-aerated volume-weighted density (ρ_(vol-weighted)) for the fluid transfer by dividing the accumulated volume-density product ({dot over (ν)}ρ_(accum)) by the accumulated volume flow ({dot over (ν)}_(accum)); and converting the non-aerated volume-weighted density (ρ_(vol-weighted)) to a standard density value using the non-aerated volume-weighted temperature (T_(vol-weighted)).
 43. The method of claim 34, with obtaining the viscosity (μ_(i)) comprising receiving the viscosity (μ_(i)) from an external source.
 44. The method of claim 34, with obtaining the viscosity (μ_(i)) comprising measuring the viscosity (μ_(i)) with the vibratory flowmeter. 